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ABSTRACT 

Typical  constructions  of  wavelets  depend  on  the  stability  of  the  shifts  of  an  underlying 
refinable  function.  Unfortunately,  several  desirable  properties  are  not  available  with  com¬ 
pactly  supported  orthogonal  wavelets,  e.g.  symmetry  and  piecewise  polynomial  structure. 
Presently,  multiwavelets  seem  to  offer  a  satisfactory  alternative.  The  study  of  multiwavelets 
involves  the  consideration  of  the  properties  of  several  (simultaneously)  refinable  functions. 
In  Section  2  of  this  paper,  we  characterize  stability  and  linear  independence  of  the  shifts  of 
a  finite  refinable  function  set  in  terms  of  the  refinement  mask.  Several  illustrative  examples 
are  provided. 

The  characterizations  given  in  Section  2  actually  require  that  the  refinable  functions 
be  minimal  in  some  sense.  This  notion  of  minimality  is  made  clear  in  Section  3,  where  we 
provide  sufficient  conditions  on  the  mask  to  ensure  minimality.  The  conditions  are  shown 
to  be  also  necessary  under  further  assumptions  on  the  refinement  mask.  An  example  is 
provided  illustrating  how  the  software  package  MAPLE  can  be  used  to  investigate  at  least 
the  case  of  two  simultaneously  refinable  functions. 
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stability  and  independence  of  the  shifts 
of  finitely  many  refinable  functions 

Thomcis  A.  Hogan 


1.  Introduction 


In  this  paper,  we  provide  a  characterization  of  the  stability  and  linear  independence 
of  the  shifts  of  finitely  many  compactly  supported  refinable  functions  $  C  L^(IR).  The 
characterization  is  in  terms  of  the  refinement  mask.  The  nature  of  these  results  leads  us  to 
a  question  involving  the  “length”  of  the  space  generated  by  these  functions.  Some  notes 
on  this  topic  are  also  provided. 

We  begin  with  some  notation  used  throughout  this  paper.  The  symbols  IN  and  7L-^ 
denote  the  set  of  natural  numbers  and  non-negative  integers  respectively.  We  also  use 

IT  ;=  [0  . .  2k) 

frequently,  but  with  the  understanding  that  any  reference  to  a  function  /  defined  on  the 
half-open  interval  TT  also  refers  to  the  27r-periodization  of  /, 

f{x  -f-  2y7r)  :=  f{x)  for  every  r  €  TT,  j  €  ZZ 
and  vice  versa.  For  any  set  S,  we  denote  by  ^^(S),  the  set  of  all  a  E  satisfying 


And  we  denote  by  (IR)  the  set  of  all  measurable  functions  /  :  IR  — C  satisfying 

ll/IU=:={/’ <«>. 

and  equipped  with  the  usual  inner  product 


dx. 


A  finite  set  $  C 
Cl  and  C2  such  that 


L^(]R)  is  said  to  have  stable  shifts  if  there  exist  positive  constants 


ci||a||f2  <  II  X]  S  <  C2||a||;2 
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for  any  sequence  a  €  x  2Z)  (it  is  often  said  that  $  provides  a  Riesz  basis  (for  its 
closed  linear  span)  in  this  case);  a  compactly  supported  finite  set  $  C  L^(IR)  is  said  to  have 
linearly  independent  shifts  if  the  map 

--EE  -  j) 

IS  one-to-one  on  U  . 

In  [JMl],  Jia  and  Micchelli  gave  necessary  and  sufficient  conditions  for  the  linear 
independence  of  the  shifts  of  a  finite  compactly  supported  set  $  C  i/^(IR).  They  showed 
that  $  has  linearly  independent  shifts  if  and  only  if  the  sequences 

(1.1)  ^ 

are  linearly  independent  for  every  ^  €  C.  In  [JM2],  moreover,  they  showed  that  a  finite 
compactly  supported  set  $  C  i/^(IR.)  has  stable  shifts  if  and  only  if  the  sequences  in  (1.1) 
are  hnearly  independent  for  every  ^  G  IR-  In  this  paper,  we  are  concerned  with  the  stability 
and  independence  of  the  shifts  of  refinable  functions. 

A  compactly  supported  finite  set  $  C  X/^(IR)  is  said  to  be  refinable  if  ($  is  not 
identically  zero  and)  there  exists  a  finitely  supported  sequence 

satisfying 

$  =  a(j)$(2  •  -i)  =  (  •  -j) 

je’a  \i€Z 

Equivalently,  $  is  refinable  if 

$(2a;)  =  A(e-*‘^)$(a;)  for  all  a;  €  C, 

where  A  :=  and 

A^,rp{z)  ^  X]  ^ 

je'z 

Unless  stated  otherwise,  whenever  we  refer  to  a  refinable  $  with  (refinement)  mask 
A,  we  mean  that 

A  =  ( 

is  a  matrix  of  trigonometric  polynomials;  that  1  is  an  eigenvalue  of  A(0)  and  a  right 
1-eigenvector  v  is  specified,  i.e., 

A(0)i;  =  V  7^  0; 
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and  that  $  C  L^(IR)  is  a  compactly  supported  solution  to  the  refinement  equation 

(1.2)  $(2-)  =  A^-,  $(0)  =  u. 

The  notions  of  stability  and  linear  independence  of  the  shifts  of  refinable  functions 
are  important  in  many  areas  including  subdivision,  finite  element  analyses,  approximation 
from  shift-invariant  spaces,  and  wavelets.  And  most  recently,  they  have  proved  to  be  im¬ 
portant  in  the  study  of  what  are  now  being  referred  to  as  multiwavelets.  Multiwavelets 
have  been  studied  by  Donovan,  Geronimo,  Hardin,  Kessler,  and  Massopust  beginning  with 
[HKM],  in  which  continuous  scaling  functions  $  were  provided.  Then  in  [GHM],  arbitrar¬ 
ily  smooth  scaling  functions  were  constructed  including  a  pair  with  orthonormal  shifts,  and 
these  scaling  functions  were  then  used  to  construct  compactly  supported  symmetric/anti¬ 
symmetric  orthonormal  multiwavelets  in  [DGHM].  Multiwavelets  have  also  been  con¬ 
structed  independently  by  Goodman  and  Lee.  In  [GL],  they  dealt  with  Riesz  bases  and 
dual  bases  for  the  set  5"°  and  its  orthogonal  complement  in  S^,  where  5°  was  a  space 
of  spline  functions  on  IR  with  multiple  knots  on  2Z  (see  also  [GLT]).  [SS]  provides  an 
alternative  pair  of  piecewise  linear  multiwavelets. 

As  with  wavelets,  it  seems  important  to  understand  the  properties  of  multiwavelets  in 
terms  of  the  refinement  mask  of  a  refinable  function  set  $.  Already  much  work  has  been 
done  in  this  area.  Approximation  orders  of  refinable  function  sets  have  been  characterized 
in  [HSS].  This  characterization  appeared  independently  in  [P],  which  also  featured  a 
factorization  of  the  mask.  This  factorization  has  already  been  put  to  use  in  [CDP]  to 
study  the  regularity  or  smoothness  of  refinable  function  sets.  All  of  these  results  were 
arrived  at  under  the  assumption  that  the  shifts  of  the  refinable  set  be  stable  or  independent. 
Yet  none  of  these  papers  provided  a  means  for  testing  the  mask  to  determine  stability  or 
independence. 

In  Section  2,  we  provide  a  characterisation  of  stability  and  linear  independence  in 
terms  of  the  mask.  This  however  brings  up  another  issue  which  does  not  appear  when 
dealing  with  a  single  refinable  function.  The  issue  concerns  the  notion  which  de  Boor, 
DeVore,  and  Ron  refer  to  in  [BDR]  as  length.  A  necessary  condition  for  a  set  $  to  have 
stable  shifts  is  that  #  be  a  minimally  generating  set.  This  unfortunately  is  not  a  condition 
which  is  easily  characterised  in  terms  of  the  mask.  For  this  reason.  Section  3  is  devoted  to 
notes  on  the  length  of  a  refinable  space  and  how  it  depends  on  the  mask. 


2.  Stability  and  independence 
2.1.  General 

This  section  deals  with  #  C  L^(IR)  which  satisfy  len5'($)  =  where 

lenS’($)  :=  min{#’i'  :  5('I')  =  5($)  }; 

S'(#)  :=  the  L^-closure  of  5o($);  and 
5'o($)  :=  the  set  of  all  finite  linear  combinations 
of  integer  translates  of  all  the  ^ 
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The  results  of  [BDR]  imply  that  for  any  compactly  supported  finite  set  $  C  T^(IR)  with 
lenS'($)  =  #$,  the  sequences  in  (1.1)  are  linearly  independent  for  almost  every  ^  G  C.  If 
len5($)  <  #$,  on  the  other  hand,  the  sequences  were  shown  to  be  linearly  dependent  for 
almost  every  ^  G  C.  Evidently,  if  len$  <  then  $  does  not  have  stable  shifts.  The 
main  results  of  this  section  characterize  stability  and  linear  independence  of  a  refinable  set 
$  in  terms  of  the  refinement  mask  A  under  the  assumption  that  len  5’($)  =  Then, 
we  complete  the  analysis  by  providing,  in  Section  3,  necessary  and  sufficient  conditions  for 
len5($)  =  in  terms  of  the  mask  A. 

Before  stating  our  main  results,  we  provide  some  additional  definitions.  For  a  given 
integer  m  >  2,  we  say  that  a  point  a;  G  IR  is  m-cyclic  in  TT  if  2’”uj  =  w  ^  0  in  TT;  i.e.,  if 
uj  ^  2E7r,  while  2'”a;  —  a;  G  27L'k.  Equivalently,  a;  G  IR  is  m-cyclic  in  TT  if  and  only  if 

2^7r 


for  some  fj.  G  2Z\(2’”  -  l)7L.  It  follows  that  G  C  is  a  primitive  odd  root  of  unity  if 
and  only  if  cj  is  cyclic,  i.e.  is  m-cyclic  in  TT  for  some  m.  We  point  out  that  if  u)  is  m-cyclic 
in  TT,  then  so  is  2^u}  for  any  £  G  2Z-|..  Also,  if  u)  is  cyclic,  then  a;  -f  tt  is  acyclic,  i.e.  is  not 
m-cyclic  in  TT  for  any  integer  m.  Finally,  it  was  shown  in  [H:  pp«  10—11]  that,  if  a?  is 
m-cyclic  in  TT  for  some  integer  m  >  2,  then  for  any  j  G  K, 

(2.1)  u;  -I-  2jTr  =  2’””a>  -|-  2’”"“^i/7r 

for  some  n  G  IN,  k  G  {0,  1,  . . . ,  m  —  1},  and  i/  G  7L\27L  (see  also  [JW;  p.  1123]). 

As  in  [H],  we  say  that  a  27r- periodic  function  /  :  C  — >  C  has  a  7r-periodic  zero 
in  IR  (resp.  C)  if  there  exists  2;  G  IR  (resp.  C)  such  that  f{z)  =  /(z  +  n)  =  0;  and  a 
contaminating  zero  if  there  exists  an  integer  m  >  2  and  to  which  is  m-cyclic  in  TT,  such 
that 

f{2'^io  -k  tt)  =  0  for  all  fc  G  {0,  1,  . . . ,  m  —  1}. 

We  also  define,  for  any  k,n  ^7L, 

n  A{2^-):=A{2"-^-)A{2^-^-)---A{2’‘+^-), 

n>e>k 

with  the  understanding  that  the  (empty)  product  represents  the  $  x  $  identity  matrix 
when  n  —  1  <  A:  -k  1. 

Lastly,  for  A  G  C^,  we  define 


A$  :=  ^  Xtjttj)  and 
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2.2.  Statement  of  main  results 


(2.2)Theorem.  Suppose  $  is  a  solution  to  the  refinement  equation  (1.2)  and  len5($)  = 
Then^  has  stable  (resp.  linearly  independent)  shifts  if  and  only  if  for  every  X  €  C^\0, 

(i)  if  A<&(0)  =  0,  then  there  exists  n  €  71+  so  that  AA”(0)A(7r)  ^  0; 

(ii)  if  AA(aj)  =  0  for  some  a;  6  H  (resp.  u>  E  C),  then  XA(u:  +  tt)  ^  0;  and 

(Hi)  for  any  integer  m  >2  and  any  u;  €  IR  which  is  m-cyclic  in  TT,  there  exists  n  €  IN  and 
k  €  {0,  1,  . . . ,  m  —  1}  so  that 


XAmnA^)M7^^  +  tt)  /  0. 

As  stated  in  Section  2.1,  if  #  has  stable  shifts  then  len5'($)  =  This  provides  the 
following  necessary  conditions. 

(2.3)  Corollary.  Suppose  $  is  a  solution  to  the  refinement  equation  (1.2).  If  #  has  stable 
(resp.  linearly  independent)  shifts,  then  for  every  A  €  C^\0, 

(i)  if  A$(0)  =  0,  then  there  exists  n  E  7Z+  so  that  AA"(0)A(7r)  ^  0; 

(ii)  if  Ayl(w)  =  0  for  some  a;  G  H  (resp.  u  E  C),  then  XA{io  +  tt)  ^  0;  and 

(Hi)  for  any  integer  m  >2  and  anyio  E  IR  which  is  m-cyclic  in  TT,  there  exists  n  €  IN  and 
k  E  {0,  1,  . . . ,  m  —  1}  so  that 


XAmn,kA)M’^^^  +  tt)  7^  0. 


It  should  be  noted  that,  in  both  of  the  conditions  (i)  and  (iii)  of  the  above  theorem 
and  corollary,  n  can  be  bounded.  Specifically,  one  need  only  verify  condition  (i)  for  n  = 
0,  1,  . . . ,  —  1;  and  one  need  only  verify  condition  (iii)  forn  =  1,  2,  . . . ,  As  regards 

condition  (i),  this  is  obvious  since  A  is  a  #  x  $  matrix.  To  see  that  n  can  also  be  bounded 
in  condition  (iii),  we  introduce,  for  given  integer  m  >2  and  iv  which  is  m-cyclic  in  TT,  the 
matrix 

M(u}')  :=  Anm+m,nm—lA). 


Since 


nnm  _  om 

2nm-ke^  _  2^+^u  =  2^— - -uj)E  2E7r 

2”*  —  1 

for  n  G  IN  and  =  0,  1,  •  •  • ,  m  —  1,  M{(jo)  is  independent  of  n  G  IN.  Moreover, 


AmnA'^)  =  ^(w)Am,fc(i*^), 


which  are  spanned  already  by  the  cases  n  =  1,  2,  •  •  • ,  for  any  particular  k.  Note  also 
that  if  det  A  is  not  identically  zero,  then  m  could  also  be  bounded  in  terms  of  the  number 
(and  locations)  of  zeros  of  det  A. 

We  also  provide  the  following  sufficient  conditions. 


5 


(2.4)  Corollary.  Suppose  $  is  a  solution  to  the  refinement  equation  (1-2)  and  lenS'($)  — 
If  A  satisfies  the  following  three  conditions,  then  the  shifts  of  $  are  stable  (resp. 

linearly  independent): 

(i)  the  matrix  ($(0)  A['n)  )  is  full  rank; 

(ii)  det  A  has  no  ir-periodic  zeros  in  IR  (resp.  C );  and 
(Hi)  det  A  has  no  contaminating  zeros. 

At  first,  the  condition  (iii)  appearing  in  (2,2)Theorem  and  (2.3) Corollary  may  seem 
to  be  nearly  inaccessible.  However,  if  one  first  tries  to  apply  (2.4)Corollary,  and  finds  that 
the  mask  does  not  satisfy  the  easily  checked  condition  (2. 4. iii),  then  one  should  begin  to 
analyze  the  (left- )nullsp aces  of  the  matrix  A  at  the  contaminating  zero  of  det  A.  Verifying 
condition  (2.2.iii)  at  these  points,  is  then  straight-forward. 

2.3.  Lemmata 

For  a  single  function  (f)  €  L^(IR),  we  refer  to  S{(j))  :=  S{{(f>})  as  the  principal  shift- 
invariant,  or  PSI  space,  generated  by  (j).  This  terminology  is  used  in  the  first  of  the 
following  results  both  of  which  can  be  found  in  [BDR]. 

(2.5) Result.  For  any  compactly  supported  Unite  set  $  C  I/^(IR),  S{^)  is  the  orthogonal 
sum  of  finitely  many  PSI  spaces,  each  generated  by  a  compactly  supported  function  having 
linearly  independent  shifts. 

(2.6) Result.  If  the  compactly  supported  finite  set  $  C  J[/^(IR)  has  linearly  independent 
shifts  and  xj)  G  S{^)  is  compactly  supported,  then 

0  -jWU) 

for  some  finitely  supported  a  € 

(2.5)Result  guarantees,  for  any  compactly  supported  finite  set  $  C  Z/^(1R),  the  exis¬ 
tence  of  a  compactly  supported  set  ^  C.  (IR)  having  linearly  independent  shifts  for  which 
5($)  =  5(^).  Then  (2.6)Result  implies  that  each  $  is  a  finite  linear  combination  of 
the  shifts  of  the  elements  of  In  the  Fourier  domain,  this  implies  the  existence  of  a  $  x 
matrix  U  of  trigonometric  polynomials  satisfying 

$  =  U^. 

If  len  S{^)  =  then  U  is  a  square  matrix.  Moreover,  since  ^  has  linearly  independent 
shifts,  the  sequences  (1.1)  are  linearly  dependent  at  ^  if  and  only  if  U{()  is  singular.  More 
specifically,  for  A  €  C^, 

(2.7)  (a$(^ -f- 2y7r)y^^  =  0  A17(^)  =  0. 

Assuming  that  len  S'($)  =  the  sequences  (1.1)  are  linearly  independent  for  almost 
every  ^  G  C.  This  implies  that  det  U  is  a  non-trivial  trigonometric  polynomial,  which  in 
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turn  implies  that  the  matrix  B,  of  rational  trigonometric  polynomials,  is  well-defined  by 
the  relation 

(2.8)  U{2-)B  :=  AU. 

This  definition  leads  to  the  relation 

i7(2-)$(2-)  =  $(2-)  =  =  AU^  =  U{2-)B^, 


which  implies  that 

(2.9)  $(2-)  =  5$ 

almost  everywhere.  The  characterization  of  S{'^)  given  in  [BDR]  then  implies  that 
'^{•/2)  €  S'(^^).  Since  the  shifts  of  axe  linearly  independent,  (2.6)Result  implies  that  the 
entries  of  the  matrix  B  axe  trigonometric  polynomials.  Now,  an  immediate  consequence 
of  (2.8)  is  that 

(2.10)  >l„,fc(a;)f/(2^+ia;)  =  U{2^uj)Bn,kH 
for  any  oj  G  C  and  integers  n  >  k  -j- 1. 

2.4.  Proof  of  main  results 

Proof  of  (2.2)Theorem.  To  see  that  condition  (i)  is  necessary  for  $  to  have 
independent  shifts,  note  that  every  element  of  22Z7r\0  has  the  form  2”‘*'^(2A;  -|-  l)7r  for  some 
n  G  ZZ^.  and  fc  G  ZZ,  and  that 

$(2"+^(2k  +  l)7r)  =  A"(0)A(7r)$((2k  +  l)7r). 

It  is  then  clear  that  if  A  does  not  satisfy  condition  (i),  then  (  A$(2j7r)  j  =  0.  Also,  if 
A  does  not  satisfy  condition  (2.2.ii),  then 

A$(2u;  -f  Ajir)  =  AA(a;)$(u;  -f-  2j7r)  =  0  and 
X^{2uj  -f  Ztt  -I-  4j7r)  =  \A{(jj  -1-  7r)$(a;  -|-  tt  -f  2j7r)  =  0  for  all  j  G  ZZ. 

So  the  sequences  (1.1)  are  linearly  dependent  at  ^  =  2a;.  Of  course  a;  G  IR  =>  ^  G  IR. 
Now,  if  A  does  not  satisfy  condition  (iii),  we  show  that 

A$  (a;  -f  2j7r)  is  zero  for  all  j  G  TL. 

Suppose  a;  is  m-cyclic  in  IT  for  some  integer  my  2  and  suppose  j  G  Z  is  given.  Then 
it  follows  from  equation  (2.1)  that 

A#(a;  2jn)  =  A$  (2’""a;  -f-  2'"”"*^i/7r) 

=  A.4„j„,jt(w)$(2'^+*w  -I-  2utv) 

=  +  I"’r)$(2''  +  i/tt). 
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for  some  n  €  IN,  fc  €  {1,  2,  dots,  m-l},ue  7L\27Zj.  So  A$(u;  +  2j7r)  is  zero  if  condition 
(iii)  is  not  satisfied. 

This  completes  the  proof  of  necessity.  To  prove  sufficiency,  suppose  that  the  shifts 
of  $  axe  dependent.  Suppose,  moreover,  that  A  satisfies  conditions  (i)  and  (ii).  We  show 
that  then  (iii)  is  violated.  ^ 

Assuming  that  the  shifts  of  $  are  dependent,  (2.7)  implies  the  existence  of  A  G  C  \0 
and  G  TT  +  ilR  so  that  Ai7(i9)  =  0.  We  will  show  that  this  implies  that  d  is  cyclic. 

Since  XU{'d)  =  0,  equation  (2.8)  implies  that 

and 

A/l  (5  +  ’')  (f  +  ’')  = 

Condition  (ii)  then  implies  that  U  is  singular  at  either  ^  or  ^  +  tt  (if  the  shifts  of  $  are 
in  fact  not  stable,  then  G  IR).  Since  det  17  is  a  trigonometric  polynomial,  it  has  only 
finitely  many  zeros  in  TT  +  zIR,  and  the  arguments  used  in  the  proof  of  [JW^:  Lemma  1] 
imply  that  -  t?  G  22Z7r  for  some  m  >  2.  We  must  still  show  that  t?  ^  2Z;7r. 

Suppose  \U{0)  =  0,  then  A$(0)  =  0  and  (i)  implies  that  AA”(0)A(7r)  is  not  zero  for 
some  n  G  E+.  Evidently, 

AA"(0)A(7r)t/(7r)  =  AC/(0)S”(0)B(7r)  =  0. 


But  this  is  a  contradiction,  since  tt  is  acyclic,  so  U{it)  is  invertible. 

Accepting  the  existence  of  the  integer  m  >  2  and  the  m-cyclic  d  for  which  \U{‘d)  is 
zero,  we  proceed  to  prove  that  this  violates  condition  (iii). 

Since 

omn  _  1 

-  u;  =  — — —  (2”*a;  -  w)  G  2E7r, 

^  X 

\U{uj)  =  0  implies  that  XU  (2’"”u;)  =  0.  By  (2.10)  and  the  27r  periodicity  of  U , 


XAmn,k{^)A(2^ U!  +  Tt)U(2^U}  +  Tt)  —  XAmn,ki^)U {2  (jj)B(2  CJ  +  Tt) 

=  Af7(2”*"a;)5,„„,jt(u;)J5(2*'a;  +  tt) 

=  0. 


Since,  2*^ a;  +  tt  is  acyclic,  U {2^u}  +  tt)  is  invertible,  so 

A^,„„,fc(cu)A(2'=c^  +  7r)  =  0 


for  every  n  G  IN,  A:  G  {0,  1,  •  •  • ,  m  —  1}. 
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Proof  of  (2.4) Corollary.  We  first  point  out  that  conditions  (i)  and  (ii)  clearly 
imply  conditions  (2.2.i)  and  (2.2.ii).  Note,  then,  that  it  is  sufficient  to  prove  the  corollary  as 
it  relates  to  stability.  For,  suppose  that  det  A  has  no  7r-periodic  zeros  in  C,  then  it  certainly 
has  no  7r-periodic  zeros  in  IR.  If,  moreover,  condition  (i)  is  satisfied  then  (2.2)Theorem 
implies  that  the  sequences  (1.1)  are  independent  for  all  ^  €  C\IR.  If  condition  (iii)  is  also 
satisfied,  and  we  accept  the  stability  statement  of  the  corollary,  then  the  shifts  of  $  must 
in  fact  be  independent. 

So,  suppose  that  (i)  and  (ii)  are  satisfied,  and  that  the  shifts  of  $  are  not  stable.  We 
will  show  that  condition  (iii)  is  violated. 

The  proof  of  (2.2)Theorem  implies  that  the  matrix  U  of  (2.7)  has  at  least  one  singu¬ 
larity,  and  that  all  such  singularities  are  at  points  which  are  cyclic.  Now,  we  show  that  for 
any  u)  which  is  cyclic, 

,  s  ,  X  .X  f  det  U  (co)  =  0  and 

{detA(a;  +  7r)  =  0. 

First,  equation  (2.8)  implies  that 

A  (cj  A  tt)  U  (u>  -h  tt)  —  U  (2io)  B  (lo  +  7r)  =0. 

And,  since  w  -|-  tt  is  acyclic  so  C/  (w  -|-  tt)  is  invertible,  this  implies  that  det  A  (u;  4-  tt)  =  0. 
Condition  (ii)  then  implies  that 

(2.12)  detA(a))^0. 

Next,  the  definition  of  B  also  implies  that 

A(u;)[/(u;)  =  l7(2co)B(u;)  =  0 

which  together  with  (2.12)  implies  that  det  17  (w)  =  0,  proving  (2.11). 

Now,  if  U  has  a  singularity  at  a  point  a;  which  is  m-cyclic  in  TT,  then  det  U(2'^u)  = 
det  U{u)  =  0.  Applying  (2.11)  repeatedly,  shows  that  det  A  has  a  contaminating  zero.  □ 


2.5.  Examples 


(2.13)Example.  This  Erst  example  will  demonstrate  the  necessity  of  the  hypothesis 
len5(#)  =  #$. 

Let  u  :=  ( 1  0  )^  and  A(u;)  ;=  A(e“*‘^)  where 


/  z^-z+2 


Then 


4(0)  =  I(1)=(J  J^) 
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satisfies  A(0)t;  =  v.  Now,  dety4(z)  =  (z  +  l)/2,  which  means  that  A  is  invertible  except 
at  TT.  Therefore,  A  satisfies  the  conditions  (ii)  and  (iii)  of  (2.4) Corollary.  Moreover, 

[$(0)  i4(7r)]  =  [u  A{-1)]  =  (J^  _2  o) 

is  invertible.  The  mask  A  satisfies  all  three  conditions  of  (2.4) Corollary.  However,  there  is 
a  solution  $  C  X/^(1R)  which  satisfies  the  refinement  equation  (1.2)  for  this  A,  which  does 
not  have  independent  (or  even  stable)  shifts. 

Let 


Then 


$  ;  = 


$(u;)  = 


2^[0..2) 


—  t2w 

—  tu; 


satisfies  the  refinement  equation  (1.2).  Moreover, 


$((2i  +  l)7r)  = 


0 

4 

(2j+l)7r 


for  any  j  €  ZZ, 


demonstrating  that  the  sequences  (1.1)  are  linearly  dependent  for  ^  —  tt.  This  implies  that 
#  does  not  have  stable  shifts.  In  this  case,  5($)  =  5(Xjo  j))-  So  lenS'($)  <  #$.  □ 

(2.14)Example.  For  this  example,  let  r  :=  ( 1  1  )^  and  A(u)  ;=  A{e  where 


A{z)  := 


1 

2 


If,  for  some  c<2  €  C,  AA(c(;)  =  0  then  A  is  a  scalar  multiple  of  (e  1).  So  if 

A>l(u;)  =  \A{uj  +  tt)  =  0,  then  A  =  0.  This  implies  that  A  satisfies  condition  (2.3.ii). 

Moreover, 

[8(0)  A(7r)]  =  (^J  ^(^2  1/2^) 

is  full-rank,  so  A  satisfies  condition  (2.3.i). 

However,  the  shifts  of  $  are  not  stable,  since  A  does  not  satisfy  condition  (2.3.iii).  To 
see  this,  note  that  ^  is  2-cyclic  in  IT.  Defining  Zq  :=  we  have 
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We  also  point  out  that 


^2n,it 


_(  M^-'^A{^)A{^)  if  A:  =  0,  and 


where 


M:=A 


27r\  _  1  /  1  Zq 

TJ"2Uo""  1 


if  fc  =  1, 


Now,  setting  A  :=  (  zq  “1 )? 


A  2''— +  7r  =0 


Moreover,  AM  =  A.  So 


for  A:  =  0,  1  and  any  n  €  IN.  This  violates  condition  (2.3.iii),  so  the  shifts  of  $  are  not 
stable. 

For  completeness,  we  point  out  that 

,  2  /Xrn  3,\ 


3  \^(i..2)/ 

is  a  solution  to  the  refinement  equation  (1.2)  and  indeed  its  shifts  are  not  stable. 


The  above  example,  in  which  det  A  is  identically  zero,  should  not  lead  one  to  believe 
that  this  property  implies  the  instability  of  the  shifts  of  $.  A  counter-example  to  such  a 
conjecture  is  provided  by 


A{z)  := 


for  which 


satisfying  the  refinement  equation  (1.2),  has  linearly  independent  shifts. 

(2.15)Example.  In  this  example,  we  use  the  results  of  this  section  to  investigate  the 
rehnable  functions  constructed  in  [DGHM].  These  were  actually  shown, in  [DGHM],  to 
have  orthonormal  shifts. 

Let  u  :=  (  \/2  1  )^  and  A(io)  :=  A(e~’^)  where 


A(z) 


s/2  f  Qs/2A  6v^z 


20  V  -1  -f  9z  -1-  9z2  -  z®  -Zs/2  -f  10\^z  -  3^^z2  J  ' 
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Then  $  is  defined  as  the  solution  to  the  refinement  equation  (1-2).  It  follows  that 


det  A{z)  =  -^(2  +  if. 


So  conditions  (ii)  and  (iii)  of  (2.4) Corollary  are  satisfied.  Moreover,  the  matrix 


[$(0)  A(7r)] 


0  4\^/5 
1  0  -8/5 


is  full  rank.  So  (2.4)Corollary  implies  that  the  shifts  are  linearly  independent  (as  long  as 
len5'($)  =  2,  which  we  will  show  in  (3.13)Example). 


3.  Notes  on  length 
3.1.  Background 

A  shift-invariant  subspace  of  T^(IR)  is  a  closed  linear  subspace  S  of  L^(IR)  which 
is  also  closed  under  shifts,  i.e.,  which,  for  each  ^  G  5,  also  contains  (/){■  —  j)  for  every 
integer  j  €  7Z.  A  finitely  generated  shift-invariant,  or  FSI  space,  is  then  a  space  of 

the  form  _ _ _ _ 

5($)  :=  span{  <;^(-  —  j)  :  </>  G  $,  i  €  2Z  } 

for  some  finite  set  $  C  T^(IR).  The  cardinality  of  a  smallest  generating  set  for  an  FSI 
space  S  is  called  the  length  of  S.  We  write 


len5:=min{#$  :  5  =  5($)  }. 


This  terminology  and  notation  is  from  [BDR]. 

In  Section  2,  stability  and  linear  independence  of  the  shifts  of  a  finite  set  of  mutually 
refinable  functions  $  were  characterized  in  terms  of  the  mask  under  the  assumption  that 
len5($)  =  #$,  i.e.,  that  the  FSI  space  generated  by  $  could  not  be  generated  by  any 
set  consisting  of  fewer  elements.  However,  no  means  for  verifying  this  assumption  by 
investigation  of  the  mask  were  provided.  In  this  section  we  give  a  characterization  of  this 
assumption  in  terms  of  the  mask. 

In  [BDR],  finitely  generated  shift-invariant  subspaces  of  ^^(IR)  were  characterized  as 
follows.  For  any  finite  $  C  I/^(IR)  and  any  /  G  L^(IR),  /  G  ■5'($)  if  and  only  if 

(3.1)  /  = 

for  some  27r-periodic  functions  r^.  They  say  that  $  is  a  basis  for  iS'(^)  if  the  functions 
in  (3.1)  are  uniquely  determined  by  /.  In  this  case,  the  functions  were  shown  to  be 
rational  trigonometric  polynomials.  They  also  prove  that  every  finite  set  $  C  I/^(IR)  of 
compactly  supported  functions  contains  a  basis  for  <S'($).  It  follows  that  for  any  finite  set 
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$  C  L^(IR)  of  compactly  supported  functions,  there  exists  C  $  and  a  $  x  ’I’  matrix  U 
of  rational  polynomials,  so  that  rank  U  =  and 

(3.2)  $(u;)  =  l!7(e“*‘^)^(u;)  for  almost  every  a;  €  IR. 

If  $  is  a  basis  for  S'($),  then  =  len  5'($).  It  follows  that  rank  U  =  #'1'  =  len  5($) 
in  (3.2).  In  fact,  for  compactly  supported  $,  $  is  a  basis  for  5'($)  if  and  only  if  len  $  = 

The  necessity  of  the  conditions  provided  in  this  section  will  not  require  that  $  C 
L^(IR).  The  background  we  provide  here  dealing  with  this  more  general  setting  is  primarily 
from  [HC]  and  [CDP].  As  with  $  C  I/^(1R.),  we  say  that  the  compactly  supported  finite 
set  #  C  X>'(IR)  is  refinable  with  mask  A  if  A  is  a  $  x  $  matrix  of  trigonometric  polynomials, 
w  is  a  given  right  1-eigenvector  of  A(0),  and  $  C  I>'(1R)  is  a  compactly  supported  solution 
to 

(3.3)  $(2-)  =  A$;  $(0)  =  v. 

Evidently, 

«=nM5?)«fe) 

7=1 

for  any  n  €  IN.  Now,  suppose  that  $  is  the  only  entire  function  satisfying  (3.3)  and  that 
the  limit 

(3.4) 

7=1 

converges  uniformly  on  compact  sets.  Clearly,  r(2-)  =  AF  and  r(0)  =  $(0).  Therefore,  F 
must  be  the  Fourier- Laplace  transform  of  $. 

When  proving  the  necessity  of  the  conditions  provided  in  this  section,  I  will  assume 
that  $  is  the  only  entire  function  satisfying  (3.3).  Sufficient  conditions  for  this  property 
are  provided  in  [HC]  and  [CDP].  The  following  theorem  from  [CDP]  provides  sufficient 
conditions  for  the  convergence  of  (3.4). 

(3.5) Result.  Suppose  A  is  a  square  matrix  of  trigonometric  polynomials.  Suppose  that 
V  is  a  right  1-eigenvector  of  A(0),  i.e.,  A(0)u  =  v  ^  0.  Suppose  p  <2,  where 

p  :=  /!)(A(0))  :=  max{  |A|  :  A(0)u  =  \u  for  some  u  ^  0}. 

Then  the  limit  (3.4)  converges  uniformly  on  compact  sets. 

Throughout  this  section,  we  use  11  to  denote  the  set  of  all  polynomials  of  one  vari¬ 
able  with  complex  coefficients  and  we  use  T  to  denote  the  field  of  (univariate)  rational 
polynomials  with  complex  coefficients,  i.e., 

:  /.»en,  s/0}. 
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Then,  for  any  set  S,  the  set 


{p:  S  ^  Ps} 


is  a  vector  space  (over  the  field  JF). 

We  interpret  any  finite  subset  P  of  also  as  a  matrix  P  €  and  vice  versa] 

and  We  identify  both  of  these  with  the  linear  map 

P  :  p^  :  u  Pu  :=  pUp 

pep 

We  make  no  distinction  between  these  three  interpretations  of  P .  For  such  P  and  any 
z  €.  C,  P{z)  C  C^,  equivalently  P{z)  €  equivalently  P{z)  :  denotes  P 

evaluated  at  z. 

When  <  oo,  the  subspace  C  P^  is  defined  for  P  C  P^  by 
P-^  :=  I  V  e  P^  :  {p,v)  '^PaVs  =  0  Vp  €  P 

I  s£S 


That  is,  P-*-  =  ker  P^,  where 

P^  :  P^  P^  V  P'^v  :=  ((p,u))pgp. 

Evidently,  rank  P  +  dim  P-^  =  dimP^  =  #5  whenever  #P  <  oo;  and  P-^-^  =  P  whenever 
P  is  a  subspace  of  p^ . 

We  denote  by  tr,  the  map  which  takes  any  function  p  defined  on  C  and  maps  it  to  the 
function  z  That  is,  for  any  function  p  with  domain  C,  crp  also  has  domain  C  and 

is  defined  by  the  rule 

(<7p)(z)  :=  p{z‘^)  for  all  z  €  C. 

Moreover,  if  P  is  any  set  of  functions  on  C,  then  aP  is  defined  by 

crP  :=  {  crp  :  p  €  P  }. 

Evidently,  crP  C  11  whenever  P  C  H  and  cr(AP)  =  (crA)((TP)  whenever  AB  is  well-defined. 

Lastly,  for  any  matrix  A  of  rational  trigonometric  polynomials,  I  denote  by  A  the 
matrix  of  rational  polynomials  defined  by 

A(e-‘‘^)  =  A{uj). 

With  these  definitions  in  place,  we  begin  our  analysis  with  a 
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(3.6)Lemma.  Suppose  <  oo.  Suppose  X  is  a  subspace  of  Then  X  has  a  basis 
U  C  n^.  Moreover,  such  a  basis  exists  with  the  property  that  rank  U(z)  =  for  every 
z  £  C;  equivalently,  det(U^U)  is  a  non-zero  constant. 

Proof.  Any  basis  of  X  consists  of  finitely  many  vectors  from  Since  #5  <  oo, 
a  basis  from  may  be  obtained  simply  by  multiplying  by  a  polynomial  which  is  a  common 
multiple  of  all  denominators  of  a  given  basis.  This  proves  the  first  assertion. 

Next  we  show  that  such  a  basis  U  C  is  everywhere  of  full  rank  if  and  only  if  the 
polynomial  p  :=  det{U'^U)  is  a  non-zero  constant.  Clearly,  p(z)  7^  0  if  and  only  if  U{z)  is 
of  full  rank.  Moreover,  the  polynomial  p  is  a  non-zero  constant  polynomial  if  and  only  if 
it  never  vanishes,  hence  if  and  only  if  U  is  everywhere  of  full  rank. 

Now,  suppose  U  is  a  bcisis  for  X.  Then  we  can  complete  U  to  a  basis  V  =:  [U  U'] 
for  Clearly,  if  V  is  everywhere  of  full  rank,  then  C/  is  as  well. 

Suppose  V  is  not  everywhere  of  full  rank.  Specifically,  suppose  3r  G  C  such  that 
det  y(r)  =  0.  Then,  det  F  €  (•  —  r)n.  Moreover,  since  V{r)  is  not  of  full  rank,  there  exists 
a  €  C  C  such  that  V{r)a  —  0.  Since  V  C  is  a  basis,  Va  7^  0  (i.e.,  Fa  is  not  the 
zero  element  of  J-^).  Therefore,  Va  £{•  —  r)n'^\0. 

If  a„  7^  0  for  some  v  £  U' ,  then  we  can  replace  v  in  U'  with  the  polynomial  vector 
V a/ {■  —r).  This  obviously  has  no  effect  on  the  basis  U  for  X,  while  it  reduces  the  degree 
of  det  V  by  one  (indeed,  det  V  has  been  multiplied  by  a„/(-  —  r)),  since  det  F  7^  0. 

If  on  the  other  hand,  a\ij,  =0,  then  there  is  some  u  €  17  for  which  a„  7^  0.  Again, 
replacing  u  in  17  with  Fa/(-  —  r),  reduces  the  degree  of  det  F  by  one.  Moreover,  although 
we  have  altered  the  basis  U,  the  resulting  set  is  still  a  basis  for  X,  since  a„  7^  0,  while 
o|  jji  —  0.  (Zl 


3.2.  Length  of  an  FSI  space 

In  this  section  we  provide  conditions  on  the  mask  A  which  are  sufficient  to  conclude 
that  len#  =  under  the  assumption  that  #  C  I'^(IR).  These  conditions  are  shown 
to  be  also  necessary  if  the  mask  satisfies  the  hypotheses  of  (3.5)Result.  The  assumption 
$  C  T^(IR)  is  not  required  to  conclude  the  necessity  of  the  conditions. 

In  the  satement  of  the  following  theorem,  we  assume  that  A  is  a  square  matrix  of 
trigonometric  polynomials;  that  1  is  an  eigenvalue  of  A(0);  and  that  u  ^  0  is  a  right  1- 
eigenvector  of  A(0).  We  assume  moreover,  that  $  C  T^(IR)  is  a  (compactly  supported) 
solution  to  the  refinement  equation  (3.3). 

(3.7)Theorem.  Iff  len5'($)  <  then  there  exists  P  C  with  =  rankP  = 
—  £,  and  such  that 

(i)  P^(1)$(0)  =  0,  while  rankP(l)  =  ^P;  and 

(ii)  (<tP^)AP-L  =  0. 

Proof.  Suppose  £  :=  len5'($)  <  Then  there  exists  ^  C  as  in  (3.2),  with 
^'^  =  £,  and  such  that 

$  =  P(e“'  )$ 

for  some  U  £  .  Moreover,  rank  17  =  £. 
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Let  P  C  be  a  basis  for  Then  =  rankP  —  —  i-  And  since  P^U  —  0, 

P^(1)$(0)  =  P'^(l)[/(1)$(0)  =  0. 

Moreover,  by  (3.6)Lemma,  we  may  assume  that  P^{1)  is  of  full  rank. 

Also, 

P^(e-‘2)AC7(e-*)$  =  P'^{e-'^)A  =  P^(e-‘2)$(2-)  =  iP^U){e-'^)^{2-)  =  0, 

again  since  P'^U  =  0.  Now,  'i'  is  compactly  supported  and  len5(’i')  =  so  ’i'  is  a 

basis  for  S'(^).  Therefore  P’^{e-'‘^  )AU{e-'  )  =  0.  Equivalently,  {aP^ AU  =  0.  Since  U 
is  of  full  rank,  it  is  a  basis  for  i/-'"'-  =  P'*",  so 

{aPfAP-^  =  0. 

□ 

To  deal  with  necessity,  we  drop  the  assumption  that  $  C  (IR) .  The  theorem  below 
is  valid  for  general  compactly  supported  distributions.  We  will  not  specify  precisely  what 
is  meant  by  5'($)  in  this  case.  However,  the  proof  of  the  theorem  below  provides,  under 
certain  assumptions  on  the  mask  A,  a  set  ^  satisfying 

(i) 

(ii)  ^  is  a  finite  linear  combination  of  the  shifts  of  #;  and 

(iii)  $  is  a  finite  linear  combination  of  the  shifts  of  'J'. 

Regardless  of  the  definitions  of  5'($)  and  len5,  it  is  reasonable  to  say  that  points  (i),  (ii), 

and  (iii)  imply  that  lenS'(#)  <  <  #$• 

In  the  statement  of  the  following  theorem,  we  assume  that  A  satisfies  the  hypotheses 
of  (3.5)Result.  That  is,  we  assume  that  1  is  an  eigenvalue  of  A(0);  that  u  7^  0  is  a  right 
1-eigenvector;  and  that  the  spectral  radius  /9(A(0))  <  2.  We  assume  moreover  that  $  is 
the  only  entire  function  satisfying  (3.3). 

(3.8)Theorem.  If  there  exists  P  C  H*  with  rankP  =  ^P  and  satisfying  conditions  (i) 
and  (ii)  of  (3.7)Theorem,  then  leniS'($)  <  —  #P. 

Proof.  We  will  provide  with  -  #P,  so  that  the  elements  of  $  are 

finite  linear  combinations  of  the  shifts  of  the  elements  of  and  vice  versa. 

By  (3.6)Lemma,  there  exists  a  basis  C/  C  H^  for  P-*-  with  the  properties  that  U{z)  is 
full  rank  everywhere  and  det{U^U)  is  a  constant.  This  second  property  implies  (say  by 
Cramer’s  rule)  that  ,  and  hence  the  left  inverse 

u-^  :=  {u'^uy^u'^, 

of  U,  also  consist  of  polynomials.  Since  rankP  =  #P,  ^17  =  —  #P- 

Now,  condition  (ii)  implies  that  AP"*"  C  (tP"*".  And  since  U  is  a  basis  for  P"*",  this 

implies  the  existence  of  P  €  for  which 

AU  =  {aU)B. 
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We  claim  that  B  also  satisfies  the  hypotheses  of  (3.5)Result. 

Evidently,  B  has  trigonometric  polynomial  entries,  since  B  =  (<tU  ^)AU .  Also, 
$(0)  €  fiy  condition  (i),  so  $(0)  =  ?7(l)v  for  some  v  €  C^.  For  this  u, 

U{l)B{l)v  =  A{l)U{l)v  =  A(1)8(0)  =  $(0)  =  U{l)v. 

And  since  17(1)  is  of  full  rank,  it  follows  that 

B{0)v  =  V. 

Now  suppose  5(0)u  =  Xu  for  some  u  ^  0.  Then  U{l)u  ^  0,  while 

A{l)U{l)u  =  U{l)B{l)u  =  \U{l)u. 

That  is,  every  eigenvalue  of  5(0)  is  also  an  eigenvalue  of  A(0).  Therefore, 

p(5(0))  <  p(A(0))  <  2. 


By  (3.5)Result, 

i=i 

converges  uniformly  on  compact  sets.  Moreover,  F  is  an  entire  function  satisfying  $(0)  = 
17(l)r(0),  as  well  as 


Cl(e-*2)r(2-)  =  Uie-'^)BT  =  AU{e-^')r. 

Since  we  assumed  that  #  was  the  only  entire  function  satisfying  (3.3),  we  must  have 

(3.9)  $  =  tl(e--)r. 

Consequently,  F  =  17“^(e“‘')$.  Since  U~^  consists  of  polynomials,  this  implies  that  F  is 
the  Fourier-Laplace  transform  of  some  compactly  supported  set  ’I'  C  X>^(1R)  which  is  a  finite 
linear  combination  of  the  shifts  of  ^  Then,  (3.9)  implies  that  the  elements  of  $  are  also 
finite  linear  combinations  of  the  shifts  of  the  elements  of  Since  , 

the  proof  is  complete. 


t 
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3.3.  Two  and  three  functions 


Suppose  ##  =  2,  then  the  results  of  Section  3.2  take  the  following  simpler  form,  in 
which  we  also  have  a  bound  for  degP.  As  we  are  assuming  that  $(0)  is  not  zero,  len5($) 
is  either  one  or  two.  For  $  C  P^(IR),  the  sufficient  conditions  then  simplify  to 

(3.10)Corollary.  Suppose  A  =:  (a;_,)  is  a  two-by-two  matrix  of  polynomials.  Suppose 
V  ^  0  is  a  right  1-eigenvector  of  A(0).  Suppose  $  C  L^(IR)  is  a  solution  to  the  refinement 
equation  (3.3).  If  lenS'($)  =  1,  then  there  exist  polynomials  pi  and  p2  satisfying 

,  ^  rj  dega2i +degp2-, 

deg  Pi  <  max{dega22, - 2 - J 

,  ^  rj  degai2 +degpi.^ 

degp2  <  max{degaii, - - - ), 


and  such  that 

(i)  pi(l)ui  +P2(l)w2  =  0,  while  one  ofpi(l)  orp2(l)  is  non-zero;  and 


00 


/  P2(z)  Q 

022(^)/ 


We  point  out  that  the  bounds  on  the  degrees  of  pi  and  p2  could  be  replaced  with  the 
simpler  bounds 


,  ^  r,  deg  ail  +  deg  021  2dega2i  +  degai2  1 

degpi  <  max  deg 022, - 2 - ’ - 3 - J 

,  ^  f,  deg  022  +  deg  ai2  deg  021  +  3  deg  012  \ 

degp2  <  max  <  degon, - - - , - - - | 

or  the  even  simpler  bounds  degpfc  <  maxjj{degaij}  which,  though  weaker  than  the  bounds 
in  the  statement  of  the  theorem,  more  clearly  display  that  the  degrees  of  pi  and  p2  can  be 
bounded  a  priori. 

Proof.  It  is  enough  to  provide  the  bounds  for  degpi  and  degp2  (the  rest  follows 
immediately  from  (3.7)Theorem).  We  may  assume  that  pi  and  p2  are  relatively  prime. 
Then  condition  (ii)  implies  that 

<Tpi  divides  a22Pi  —  a2iP2?  and 
crp2  divides  aiiP2  —  ai2Pi. 

The  bounds  on  degpi  and  degp2  follow.  □ 

The  results  of  Section  3.2  take  a  simple  form  in  the  czise  =  3  as  well.  Now  there  are 
two  possibilities,  lenS'(#)  =  1  or  len5($)  =  2.  Note  that  the  three  equations  in  condition 
(ii)  of  the  next  two  corollaries  are  inherently  redundant.  However,  since  we  don’t  know 
which  if  any  of  the  pk  might  be  zero,  we  have  no  choice.  To  handle  the  possibility  that 
len5($)  =  1,  we  offer  the 
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(3.11)Corollary.  Suppose  A  =:  (aij)  is  a  three-by-three  matrix  of  polynomials.  Suppose 
V  0  is  a  right  1-eigenvector  of  A(0).  Suppose  $  C  iy^(IR)  is  a  solution  to  the  refinement 
equation  (3.3).  IflenS(^)  —  1,  then  there  exist  polynomials  pi ,  p2  and  ps  satisfying 


1 


,  ^  f,  degai2  +  degp2  degoia  +  degps  1 

degpi  <  max  deg  an, - - - , - - - > 

j  ^  f  deg  021  +  deg  Pi  deg  023  +  degpa  \ 

degp2  <  max< - - - , deg 022, - 2 - f 

f  deg  031  +  deg  Pi  deg  032+ deg  p2  ,  \ 

degp3  <  max  < - - - , - - - ,  deg  033  > 


and  such  that 

(i)  P2(1)U3  -  P3(1)U2  =  P3(1K  -  Pi(l)u3  =  P\{1)V2  “  P2{l)v\  =  0,  while  not  all  of  pi(l), 
P2(l),  and  P3(l)  are  zero;  and 

00 


/  0  -P3(z^)  P2(z^)  \  f an{z)  a^iz)  ai3(z)\  / Pi{z)\ 

Pziz'^)  0  -Pi{z'^)  a2i(^)  022(2)  023(2)  P2(2)  =0. 

\-P2{z^)  Pl{z^)  0  /  \03i(2)  032(2)  033(2)/  \p3{z)J 

Again  we  point  out  that  the  bounds  above  imply  that  degp*  <  maxjj{deg  Oij}. 

Proof.  Assuming,  as  we  may,  that  the  polynomials  pi,  p2,  and  pz  are  relatively 
prime,  condition  (ii)  implies  that 

opi  divides  onPi  +  012P2  +  ciizpz, 
ap2  divides  a2ipi  +  a22P2  +  023JP3>  and 
apz  divides  osipi  +  az2P2  +  I33P3. 

The  bounds  given  now  follow  easily.  The  rest  follows  from  (3.7)Theorem.  □ 

To  handle  the  possibility  that  len5'($)  =  2,  we  make  the  additional  assumption  that 
det  A  is  not  identically  zero.  For  this  case,  we  offer  the 

(3.12)Corollary.  Suppose  A  =:  {aij)  is  a  three-by-three  matrix  of  polynomials,  for  which 
det  A  is  not  identically  zero.  Suppose  v  ^  0  is  a  right  1-eigenvector  of  A(0).  Suppose 
$  C  Z/^(IR)  is  a  solution  to  the  rehnement  equation  (3.3).  JflenS'($)  =  2,  then  there  exist 
polynomials  pi,  p2,  and  pz  satisfying  degpk  <  2maxj,j  dego^j  for  k  =  1,  2,  and  3,  and 
such  that 

(i)  Pi(l)ui  +P2{1)V2  +P3(1)V3  =  0,  while  one  of  pi(l),  P2(l),  orp3(l)  is  non-zero;  and 

00 


013(2)  \ 
023(2)  1 
033(2)  / 


/  0  -P3(2)  P2(2)  \ 

P3(2)  0  -Pi  (2)  =0. 

\-P2(2)  Pi  (2)  0  / 


(011(2)  Oi2(z) 
021(2)  022(2) 

031(2)  032(z) 
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Proof.  _ 

Let  Ma  •=  be  the  matrix  of  cofactors  of  A.  That  is, 

Since  det  A  is  not  identically  zero,  A  is  almost  everywhere  invertible,  and 

MJ  =  (det 

Moreover,  degm,j  <  2max,/j'{degai'j'}  for  i,j  =  1,2,3. 

By  (3. 7)Theorem,  there  exists  =:(  Pi  P2  pa  )  such  that 

aP^A  =  uP^ 

for  some  u  =:  ffg  £  P.  This  in  turn  implies  that 

gdetA^pT  ^ 

As  we  may  assume  that  pi ,  p2 ,  &n-d  pa  are  relatively  prime  and  as  the  elements  of  Ma 
polynomials,  /  must  divide  pdet  A  and,  in  turn,  aP^  divides  P'^Ma-  This  then  implies 
that  every  element  of  P  is  of  degree  less  than  max,  j{degm,j}  <  2maXi,j{dega,j}.  □ 

3.4.  Example 

(3.13) Example.  We  complete  the  analysis  of  the  rehnable  functions  from  [DGHM]  be¬ 
gun  in  (2.15)Example. 

Recall  that 

~  6\/2  +  6\/2z  16  \ 

“  20  V  -1  +  9z  +  9z^  -  z^  -3y/2  +  10y/2z  -  3y/2z^  )'  V  ^  / 

By  (3.10)Corollary,  if  len5($)  =  1,  then  there  would  be  polynomials  pi  and  p2  satisfying 

(3.14)  {Pi{z^)  1^2(2^) ) ^(^)  =0 

as  well  as  \^pi(l)  +  P2(l)  =  0.  The  bounds  provided  there  for  the  degrees  of  these 
polynomials  implies  that  degpi  <  2  and  degp2  <  1.  We  use  Maple  to  try  to  find  such  pi 
and  p2.  We  let 

pi{z)  =:  Go  +  aiz  +  02^;^  and  P2(^)  =: 

The  following  Maple  code  then  sets  up  these  two  conditions  and  attempts  to  solve  for  the 
coefficients  oto,  02,  bo,  61. 
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>  A;=sqrt(2)/20*matrix( [[6*sqrt(2)*(l+z) , 16] , 

>  [-l+9*z+9*z''2-z"3,  sqrt(2)*(-3+10*z-3*z''2)]] )  ; 

>  v:=matrix( [[sqrt (2)] , [1]] ) ; 

>  pl:=a0+al*z+a2*z‘'2;  p2:=b0+bl*z; 

>  sPT:=matrix([[subs(z=z“2,pl)  ,subs(z=z''2,p2)]] )  ; 

>  Pperp:=matrix([[p2] , [-pi]]) ; 

>  condi:=subs(z=l,evalm(sPT&*v)) [1,1] ; 

>  exprii :=evalm(sPT&*A&*Pperp) [1,1] ; 

>  condii:=coeffs (collect (exprii.z) ,z) ; 

>  solve({condi,condii},{a0,al,a2,b0,bl}) ; 

Maple’s  response  is 

{6o  =  0,  <22  “  0)  =  0,  (2i  =  0,  Oq  =  0} 

{fej  =bi,  02  =  —  V^&i,  <2o  =  “2'^’  =  0,  60  =  ~&i} 

Corresponding  to  pi  =  p2  =  0;  and  pi  =  ^{z'^  -  1),  p2  =  z  —  1,  respectively.  However, 
neither  of  these  choices  for  pi ,  p2  satisfies  the  condition  that  ( pi  (1)  p2(l) )  he  of  full  rank. 

Therefore,  len5($)  7^  1. 

It  is  worth  noting  that  the  polynomials  pi{z)  =  ^(1  +  z)  and  P2{z)  =  1  satisfy 
condition  (3.14).  Moreover  P^(l)  =  (pi(l)  P2(l))  =  ( 1)  has  full  rank.  However, 
P'^(l)v  =  3  7^  0.  Multiplying  pi  and  p2  by  z  —  1  does  not  effect  (3.14)  and  clearly  makes 
P^(l)v  =  0.  But  in  this  case,  P^{1)  =  (0  0)  is  not  of  full  rank.  So  this  example  also 

illustrates  that  all  three  conditions  are  necessary  to  conclude  that  len5($)  <  □ 
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